Detection of neutralino annihilation photons from external galaxies 

E. A. Baltz 

Department of Physics, University of California, Berkeley, CA 94720, USA 

C. Briot, P. Salati, and R. Taillct 
Laboratoire de Physique Theorique LAPTH, BPllO, F-74941 Annecy-le-Vieux Cedex, France 

J. Silk 

Department of Astrophysics, Oxford University, Oxford 0X1 3RH, United Kingdom 



Draft Version of February 1 , 2008 

We consider neutralino annihilation in dense extragalactic systems known to be dominated by 
dark matter, in particular M87 and several local dwarf spheroidal galaxies. These annihilations can 
produce energetic gamma rays which may be visible to atmospheric Cerenkov telescopes. We explore 
the supersymmetric parameter space, and compute the expected flux of gamma-rays coming from 
these objects. It is shown that some parts of the parameter space lead to a signal observable with 
the next generation of Cerenkov telescopes, provided the supersymmetric dark matter has a clumpy 
■ structure, as may be expected in a hierarchical scenario for structure formation. 
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I. INTRODUCTION 



^\ • One of the favorite candidates for the astronomical missing mass is a neutral weakly interacting particle. Such a 
I species is predicted in particular by supersymmetry, a theory that is actively tested at accelerators. It is conceivable 
0^ ■ therefore that most of the dark matter in the halo of the Milky Way is made of such particles. The mutual annihilation 
■ of these so-called neutralinos would yield, among a few other indirect signatures, a flux of high-energy gamma-rays. 
The latter has been extensively studied in the literature [1-4]. It is unfortunately spoiled by the diffuse background 
produced by the spallation of cosmic-ray protons on the interstellar gas as discussed by Chardonnet et al. [5]. The 
Q , distribution of molecular hydrogen inside the galactic ridge is not sufficiently well known to ensure an accurate 
5^ • prediction of the diffuse emission. This may obliterate a reliable interpretation of any putative gamma-ray excess in 
' terms of supersymmetric dark matter. 



. . , We examine here the possibility of observing the gamma-ray signal from extra-galactic systems that contain large 
' amounts of unseen matter. The giant elliptical galaxy M87, at the center of the Virgo Cluster, provides an excellent 
k>( I illustration. It is known to contain ten times more mass than our own Milky Way. Actually, X-ray observations of 
■ its ambient hot gas indicate that in the innermost 100 kpc, the mass reaches ~ 10^^ Mq. If neutralinos pervade this 
^ galaxy, a strong gamma-ray emission should be produced in the central region. Neutralino annihilation is associated 
to the gamma-ray flux at the Earth 

<^>T' = ^MA^7 / <ds (1) 

J los 

which, in the case of M87, is two orders of magnitude larger than for own galaxy. On the other hand, M87 is quite 
distant, about 15 Mpc away. It should therefore appear as a bright gamma-ray spot on the sky, extending for at most 
a few hundred square arcminutes. The detection of such a strong but localized source is well-suited for atmospheric 
Cerenkov telescopes (ACT) which can only monitor small portions of the sky at the same time but have very large 
effective collecting areas. 

Dwarf spheroidal (dSph) galaxies are also suspected to contain significant amounts of dark matter. They orbit round 
the Milky Way and are intermediate in size between the globular clusters and larger systems such as the Magellanic 
clouds. The typical mass of dSph's is ~ 1 — 4 x 10'' Mq while their radii reach up to ~ 1 — 2 kpc. They are also closer 
than M87. 

The gamma-ray flux (1) produced by annihilating species in M87 may be split into an astrophysical piece and a 
particle physics piece. The former term merely amounts to the integral along the line of sight of the dark matter 
density squared. It is approximately given by where M denotes the mass contained within the radius R. In 
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the case of M87, that Hne of sight integral is ^ 10 M©^ pc^'"^ whereas it is ^ 0.01 — 1 M©^ pc^^ for dwarf sphcroidals. 
This compares to the Milky Way value of ^ 0.1 M©^ pc^^ which corresponds to the irmer 100 kpc. These rough 
estimates arc improved in Sec. II where the distribution of matter inside M87 and the dSph systems is modeled 
together with the radial profile of the neutralino induced gamma-ray emission. The various backgrounds to the latter 
are discussed in Sec. Ill where a fiducial example is presented. In particular, the signal-to-noise ratio is derived as a 
function of the angular distance to the centers of the sources under scrutiny. The supersymmetric model is presented 
in Sec. IV. We then delineate in Sec. V the domain of neutralino masses and gamma-ray production cross sections 
which the next generation of atmospheric Cerenkov telescopes will explore. Both continuum and monochromatic 
channels are featured. We finally draw some conclusions in Sec. VI where we pay special attention to a possible 
clumpy structure of the neutralino distribution. 



II. THE M87 AND DSPH SYSTEMS 



The observed X-ray emission produced by clusters of galaxies such as Virgo results from the thermal bremsstrahhmg 
that takes place in the hot diffuse intracluster gas (see [6] for a complete review). The temperature profile T{r) of 
the gas can be inferred from the spectrum of the X-ray radiation. The electron number density profile ne(r) can be 
obtained from the X-ray surface brightness 

S(r) = Jne{Vr^ + s^)Lxds , (2) 

where Lx denotes the X-ray luminosity of the plasma. The properties of the hot gas in the vicinity of M87 have been 
determined by Tsai [7]. The X-ray data are well fitted by 

with no = 4.31 X 10-2 cm-3, ai = 6.63 kpc, as = 4.58x10^ kpc, ai = 0.49, as = 0.114, = 0.869 andToo = 8.35x10^ 
K. Assuming that the intracluster gas is in hydrostatic equilibrium, the total mass profile is given by 

M(r) = --^f^ + ^) , (4) 
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where nrrip is the average mass of a particle in the gas. The total mass inside M87 reaches 8.4 x 10^^ M© at 10 kpc 
and 1.4 x 10^^ M© at 100 kpc. The contribution from the stars alone was obtained from B-band surface brightness 
profiles (see [7] and references therein). The dark matter component distribution is derived once the stars and the 
gas have been removed. Mappings of the X-ray emission in the Virgo cluster are consistent with the hypothesis that 
this system and its central galaxy contain non-baryonic dark matter. The ratio of the dark matter density to the 
baryonic matter density is also borrowed from [7] 

0.15 f-^V ' . (5) 
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The integral along the line of sight of p^j^ is readily obtained. 

Dwarf spheroidal (dSph) galaxies also seem to contain large amounts of unseen material [8]. Like globular clusters, 
dSph's have low luminosities of order 10^ — lO'' L© within an ill-defined radius rt of a few kpc. Their structure also 
follows a King profile (see [9] for a review). If these systems are in virial equilibrium, a velocity dispersion cr of a few 
km/s translates into a mass 

M.s,..2.3xlO'M,(^)^(liE2) , (e) 

hence these systems have mass-to-light ratios that may reach ~ 200. We have been interested here in four galaxies 
with exceptionally high M/L-values [10] and galacto centric distances less than 100 kpc. The dSph's featured in table I 
are Carina, Draco, Ursa Minor and Sextans. We have modeled their mass distributions with a one-component King 
profile [11] for which the phase space density is given by 
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L M/L M d (Jobs c rt crth 

Name (10^ Lq) (Mq/Lq) (10^ M©) (kpc) (km/s) (kpc) (km/s) 



Carina 2.4 ± 1.0 59 ±47 14.1 ± 10.0 85 ± 5 6.8 ±1.6 

Draco 1.8 ±0.8 245 ± 155 44.1 ±24.0 72 ±3 10.2 ±1.8 

Ursa Minor 2.0 ± 0.9 95 ± 43 19.0 ± 10.5 64 ± 5 12.0 ±2.4 

Sextans 4.1 ± 1.9 107 ± 72 43.8 ± 23.6 83 ± 9 6.2 ± 0.8 



TABLE I. These four dSph galaxies liavo very large M/L-values and are within 100 kpc from the Milky Way center. The 
luminosities L, mass-to-light ratios, total masses M, galactocentric distances d and concentrations c have been borrowed from 
[10]. The velocity dispersions Uohs are those given by [13]. These systems have been modeled with King profiles as explained in 
the text. Consistency may be checked by comparing the resulting velocity dispersions Uth with the observed values aoha- Also 
indicated are the derived tidal radii rt, in kpc. 

f{E) cc exp (^^) - 1 . (7) 

The stellar or particle energy per unit mass is denoted hy E = + while (pt = (pi^t) is the potential at the tidal 
radius rt of the dSph. Beyond that boundary, stars have enough energy to escape from the gravitational field of the 
system and are captured by the tidal field of the nearby Milky Way. The mass distribution ensues from the phase 
space density (7) 

P{i') ^ ^exp(u)erf(Vw) — , (8) 

where u denotes the ratio {(pi — (p{r)} ja^ while erf is the error function. One component King models are completely 
determined once the velocity dispersion tr, the central density pc and the ratio x = {</>* — '/'(O)} jcP' are specified. Any 
combination of these quantities suffices. The concentrations of the four dSph's under scrutiny are given in table I as 
determined by [10]. They are related to the tidal and core radii through 

c = logio f , (9) 



r, 



where the core radius is defined as 

re = . (10) 

Values of the concentration of 0.5 and 1 respectively translate into the ratios % = 1.97 and 4.85. Fixing the concentra- 
tion allows for the determination of the density profile p{r)/pc as a function of the reduced radius r/vc- The second 
input to our calculations is the dSph mass as given by the central values of table I. We finally set the average dSph 
density pdSph by requiring that the proper gravitational field of these galaxies is compensated, at their boundaries, 
by the tidal field of the Milky Way taken at perigalacticon. Assuming a logarithmic galactic potential implies that 

PdSph = /(e) PG , (11) 
where the function /(e) depends on the orbital eccentricity of the satellite galaxy [12] 

1 + [(1 + ef /2el In [(1 + e) / (1 - e)] 

m = 77^2 • (12) 

(1 - e) 

The average galactic density pc is understood within the sphere centered on the Milky Way, with radius equal to the 
semimajor axis a of the dSph orbit. The corresponding mass is given by 

Mg = 1.1 X 1O^°M0 X a [kpc] . (13) 
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We set the semimajor axis a equal to the dSph galactocentric distance d. Eccentricities in the range between and 
1/2 are generally assumed. Here, their values were derived by requiring that the resulting velocity dispersions crth 
should be close to the observations - see table I. In the case of Carina and Draco, the eccentricity is e = while 
for Ursa Minor, a value of e = 1/2 is appropriate. As regards Sextans, we simply imposed the equality between the 
dSph and galactic average densities, i.e., pdSph = Pg, in order to get a velocity dispersion of 8.8 km/s, not too far 
from the measured value of 6.2 ± 0.8 km/s. Assuming relation (11) and a circular orbit would have lead to Cth = 9.9 
km/s. In modeling the dSph galaxies, we have adopted the point of view that these satellites contain large amounts 
of dark matter as implied by their high velocity dispersions. In the alternative explanation, the latter are the result 
of an ongoing tidal stripping. However, Oh et al. [13] showed that unbound but not yet dispersed systems have 
velocity dispersions fairly close to the values derived from the virial equilibrium. The large M/L- values of the four 
dSph systems at stake do not depend therefore on whether these galaxies are in equilibrium or are currently being 
torn apart. Notice finally that most of the putative neutralino-induced gamma-ray signal emitted by these dSph's is 
produced in the inner 10 arcmin as shown in the right panel of Fig. 1. It weakly depends on the actual boundaries of 
these systems. 



III. THE GAMMA-RAY SIGNAL AND ITS BACKGROUNDS 



The gamma-ray flux produced by a halo of dark matter particles may be computed as follows. First, for each set 
of supersymmetric parameters, we derived the thermally averaged annihilation cross section {av)A into some channel 
A. Then, we used the Lund Monte Carlo [14] to compute the photon spectrum per annihilation from hadronization 
or decay processes, and we integrated over some specified threshold which depends on the atmospheric Cerenkov 
telescopes in question. Finally, we have summed over all channels to find the total photon spectrum from a given 
supersymmetric model. The volume production rate of gamma-rays is simply given by n'^{av)Nj, where n^- is 
the number density of annihilating particles and is the mean number of photons above the threshold. The 
corresponding flux at the Earth depends on the integral, along the line of sight, of the dark matter density squared 



1 (av)N, 

as well as on the gamma-ray production cross section (av) and the mass m-^ of the species. 



Atmospheric Cerenkov telescopes (ACT) offer the most promising method to detect gamma-rays with energies 
above 100 GeV. When passing through the atmosphere, incoming high-energy photons emit a Cerenkov radiation 
that is detected and analyzed to infer the energy and arrival direction of the primary gamma-ray. The performances 
of ACTs arc determined by the surface of the optical reflectors used to collect the Cerenkov light, and the ability to 
reject the cosmic-ray induced background. A thorough analysis of these performances can be found in Aharonian [15]. 
We will simply model them by the effective collecting area A^g, a feature of the instrument that takes into account 
in particular the actual surface of detection together with the energy dependent detection efficiency. 

The surface brightness fi^^^ of the source as seen by an ACT depends on the collecting area of the instrument and 
on the integration time T 

^susy ^ ^---YA^sT . (15) 

We plot the surface brightness in Fig. 1 for a collecting area A^s = 1 km^ and for an entire year of effective observation. 
In the fiducial model taken here, the neutralino mass is my- = 1 TeV while the annihilation cross section has been set 
equal to {av) Nj = 10~^^ cm^ s~^. The surface brightness corresponds to the number of photons above a threshold 
of 100 GoV that are produced per square arcmin of the source. In the left panel, the giant galaxy M87 is fc^aturcd 
whereas in the right panel, the four dSph satellites discussed in the previous section, i.e., Carina, Draco, Sextans and 
Ursa Minor are presented. The signal declines away from the centers of these systems. It follows the integral along the 
line of sight of the neutralino density squared - remember that n-^ = p^/m-^. Actually, we must take into account the 
finite angular resolution of the instrument. When the telescope is pointed towards the center of the object, photons 
coming from directions within the resolution cone also reach the detector. Thus, the relevant quantity to examine is 
the integral of the surface brightness ^^^^ over a disk centered on the source, with angular radius 9 
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AT, = / 2naif^''^da . (16) 
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FIG. 1. The radial profiles of the neutralino-induced signal (solid curves) and of the various backgrounds (dotted and 
dashed lines) are plotted as a function of the angular distance to the source centers. A fiducial model with rrix = 1 TeV 
and (av)cont.N^ = 10~^^ cm"^ s~^ is taken while a threshold of 100 GeV is assumed. In the left panel, the backgrounds 
are respectively labeled as (a): electronic; (&): hadronic; (c): extragalactic; (d): M87 and (e): Milky Way gamma-ray diffuse 
emissions. The right panel shows nearby dwarf spheroidal galaxies with exceptionally high M/L-values. The neutralino-induced 
signal is flat within 10 arcmin from the center. The total background is featured by the dashed line (g). Carina is the exception 
with a slightly larger value (f) dominated by the Milky Way diffuse emission. That system lies within the galactic ridge. 



When the angular resolution of the ACT is good enough, the size of that circular region is generally specified by the 
requirement that the signal-to-noise ratio should be optimal. 

The mimbcr of ncutralino induced source photons Ng should actually be compared to the background A'bg. the 
origins of which are various. Would it be constant, this background could easily be removed from the signal. Actually, 
it follows Poisson statistics, and as such, it exhibits fluctuations of amplitude ^/Nh^■ These fluctuations have a smaller 
chance of being interpreted as a signal when the significance 




is large. There arc two main types of background events. First, an experimental background is due to hadronic and 
electronic cosmic-rays that impinge on the top of the atmosphere. The induced showers can be misinterpreted as 
gamma-ray events. Electrons make up the largest source of background insofar as their showers are of the electromag- 
netic type and cannot be disentangled from those generated by the impact of high-energy photons. The corresponding 
flux steeply decreases at high energy [16] 

/ rp \ -3.3±0.2 

$e = (6.4 X 10-2 GeV-^cm-^s-^sr-i) (yq-^) ' (18) 

and leads to the background 

^(E > E.) = (1.9 X 10' km- yr- ^cmi.,-) (^) . (19) 

In the left panel of Fig. 1, that electron induced background is featured by the dotted line (a). It is more than an order 
of magnitude larger than the neutralino-induced signal at maximum. It is also noticeably flat all over the source. 
Observations performed between 50 GeV and 2 TeV yield a hadron flux [17] 
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$had = (l.8GeV-icm-2s-isr-i) (j^) ' (20) 

Hadron-induced showers are more extended on the ground than those of the electromagnetic type. Stereoscopy is 
a powerful tool to discriminate hadrons from electrons and gamma-rays. The CAT experiment, for instance, has 

already achieved a rejection factor of one misidcntificd event over a sample of 600 showers generated by cosmic -ray 
hadrons [18]. Here, we have assumed an even better rejection factor, with only one misidentified hadron out of a 
thousand showers. This yields a hadron background 

> = (8.7 X 10> km- y.-' „cmi„-) {^f"' . (21) 

which corresponds to the dotted line (b) of Fig. 1. Once again, that background is flat over M87. Next come the 

astrophysical sources of background. To commence, Srcekumar ct al. [19] have measured an extragalactic component 
in the gamma-ray diffuse emission with the EGRET instrument on board the Compton gamma-ray observatory 

/ x -2.10±0.03 

$eg = (7.32 ± 0.34 X 10"^ MeV"^ cm'^ s"^ sr'^) ( ) . (22) 



This translates into the background (c) 

^iE>Eo) = (2.1xl02km-2yr-iarcmin-2) (^^^^) ' ' • (23) 

In our fiducial example, the neutralino-induced signal exceeds the extragalactic background within ~ 1 arcmin from 
the center of M87. Then, we have modeled the gamma-ray diffuse emission from that giant elliptical galaxy itself. 

Local cosmic rays interact with the chistcr gas to produce high energy photons that may potentially contaminate the 
signal. In the inner 10 kpc, the magnetic field of M87 is comparable to that of the Milky Way, with a magnitude ~ 1 
/iG. It falls by a factor of ten outwards in the Virgo cluster, at a distance of ~ 100 kpc. Assuming that equipartition 
of energy holds between this magnetic field and the cosmic-rays that pervade M87 - just like in our own galaxy - we 
infer a gamma-ray emissivity of 



Iu{E) = (2xlO-3^GeV-is-isr-i) (^^) 



-2.73 

(24) 



per hydrogen atom illuminated by high- energy protons. Once multiplied by the hydrogen column density across M87, 
it yields the in situ diffuse gamma-ray flux. We have assumed that hydrogen is the main constituent of the gas at the 
center of the Virgo cluster so that we have integrated the electron density, as derived in Sec. II, along the appropriate 
lines of sight. The gamma-ray emissivity as given by relation (24) is typical of the inner 10 kpc inside M87. It has 

been rescaled by the factor |l + (r/10kpc)^| to account for the outward decrease of both the magnetic energy 

and cosmic-ray flux. The M87 gamma-ray diffuse background is presented as the dashed line (d) of Fig. 1. It falls 
outwards as a result of the combined decrease of the hydrogen column density and of the cosmic-ray fltix. Finally, we 
have taken into account the Milky Way diffuse emission. We use hydrogen column densities inferred from the dust 
maps of [20]. A hydrogen column density of 1.69 x 10^'^ H cm~^ in the direction of M87 corresponds to 

> E., = (2.8 (^)"" . (25) 

and to the dashed line (e). The Milky Way component is the weakest source of background. Electrons and in a lesser 
extent hadrons - as long as rejection is efficient - dominate. In the right panel of Fig. 1, the total background is 
presented. It encompasses the various sources discussed above except a local diffuse emission. The dSph galaxies 
contain old stars and the acceleration mechanisms at work in the Milky Way are presumably less prevalent there. Even 
in the case of the brightest source, Ursa Minor, the signal is still three orders of magnitude below the background - see 
curve (g). Because Carina lies in the sky towards the galactic ridge, the hydrogen column density of the Milky Way 
is quite large, reaching a value of 8.7 x 10^^ H cm~^. The gamma-ray diffuse emission of our own galax;y dominates 
the background (f) in that specific case. 

The significance S is presented in Fig. 2 for M87 (left panel) and in the case of Ursa Minor, the best dSph source 
(right panel). Should the signal be flat, both Ns and A'bg would be proportional to the surface of the sky monitored by 
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FIG. 2. Signal-to-noisG ratio as a function of threshold and beam size. Again a fiducial model with — 1 TeV and 
(o"^')cont. Nj = 10^^'' crn"^ s^^ is taken, while a more realistic integration of 0.01 km'^ yr is now assumed. The left panel 
features the case of M87 for which the significance 5 peaks at an angular radius of ^ 1.4 arcmin. For the brightest dSph Ursa 
Minor (right panel), the optimal beam size is 7.7 arcmin. Three values of the threshold energy are presented. 



the ACT. As a result, the signal-to-noise ratio would just increase like the angular radius 9 of the observed region. In 
Fig. 2, it actually increases as the beam size opens up. Because the signal is not flat but weakens far from the source, 
the signal-to-noise ratio reaches a maximum and drops at larger radii. The ACT acceptance has been set equal to 0.01 
km^ yr. The neutralino mass is still 1 TeV whereas the gamma-ray production cross section {av)cont. = 10~^^ 
cm"* s^^ is the same for the three values of the threshold energy featured here. As the background is flat over the 
source, the signal-to-noise ratio always peaks at the same angular position. For M87, the significance is the largest 
for an angular aperture of ~ 1.4 arcmin. Because Ursa Minor is closer, the optimal beam size becomes 7.7 arcmin. As 
both the electronic and hadronic backgrounds weaken with the gamma-ray energy, the magnitude of the maximum 
increases with the threshold as is clear in both panels. For M87 and a 50 GeV threshold, the significance reaches a 
value of ~ 0.4 whereas for 250 GeV, the peak value increases up to 2. In the case of Ursa Minor, the variations of the 
significance are just the same. Should the acceptance be 1 km^ yr, these signal-to-noise ratios would be rescaled up 
by an order of magnitude. 



IV. EXPLORING THE SUSY PARAMETER SPACE 



We have explored the Minimal Supersymmetric Standard Model (MSSM). This framework has many free param- 
eters, but with reasonable assumptions the set of parameters is reduced to seven: the Higgsino mass parameter /z, 
the gaugino mass parameter M2 , the ratio of the Higgs vacuum expectation values tan /3, the mass of the CP odd 
Higgs boson toa, the scalar mass parameter mo and the trilinear soft SUSY-breaking parameters Ai, and At for third 
generation squarks. The only constraint from supergravity that we imposed is gaugino mass unification, though the 
relaxation of this constraint would not significantly alter our results. For a more detailed description of these models, 
see Refs. [21,22]. 

The lightest stable supersymmetric particle in most models is the lightest of the neutralinos, which are superpositions 
of the superpartners of the gauge and Higgs bosons. We used the one-loop corrections for the neutralino and chargino 
masses given in Ref. [23]. For the Higgs bosons we used the leading log two-loop radiative corrections, calculated 
within the effective potential approach given in Ref. [24]. We made thorough random scans of the model parameter 
space, with overall limits of the seven MSSM parameters as given in Table H. 

Each model was examined to see if it is excluded by the most recent accelerator constraints. The most important 
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TABLE II. The ranges of parameter values used in our scans of the MSSM parameters. In total, we used approximately 



116,000 models that are not excluded by current accelerator constraints. 



of these are the LEP bounds [25] on the hghtest chargino mass 



,91 GeV , |m^.-m,o|>4GeV 



\ 85 GeV , otherwise 

and on the hghtest Higgs boson mass m^o (which ranges from 72.2-88.0 GeV depending on sin(/3 — a) with a being 
the Higgs mixing angle) and the constraints from 6 — > 57 [26]. 

For each allowed model, the relic density of neutralinos O^/i^ was calculated, where Sl^ is the density in units of 

the critical density and h is the present Hubble constant in units of 100 km Mpc"-'^. We used the formalism 
of Ref. [27] for resonant annihilations, threshold effects, and finite widths of unstable particles and we included all 
two-body tree-level annihilation channels of neutralinos. We also included the so-called coannihilation processes 
between all neutralinos and charginos in the relic density calculation according to the analysis of Edsjo and Gondolo 
[21]. Present observations favor /i = 0.6 ± 0.1, and a total matter density Om = 0.3 ± 0.1, of which baryons may 
contribute 0.02 to 0.08 [28]. However, wc only required that il^h^ < 0.5. We were also interested in models where 
neutralinos are not the only component of dark matter. In models with ^^h^ < 0.025, we rescaled the relevant halo 
densities by a factor 0^^/1^/0.025 to account for the fact that a supplemental source of dark matter is required in such 
models. 



V. FLUX OF GAMMA-RAYS FROM NEUTRALINO ANNIHILATION 

We now present the results of the scans over the supersymmetric parameter space. We first made scatter plots of the 
continuum gamma-ray production rate (cri')cont. versus the neutralino mass m^, for three gamma-ray thresholds 
Eth ~ 50, 100 and 250 GeV. The results are presented in Fig. 3, along with the 3-it detection limit of the signal from 
M87, for a typical exposure of 0.01 km^ yr. This acceptance corresponds to the next generation of ACTs. The HESS 
project [29] for instance should reach a threshold of 40 GeV for a collecting area of 300 m by 300 m. Its angular 
resolution should be 0.1°. We have assumed a generous 0.1 yr integration time which would correspond to a few 
months of continuous observation. This compares to the VERITAS project which should be sensitive to the energy 
range extending from 50 GeV up to 50 TeV. A collecting area of 10,000 m^ is expected at 100 GeV, increasing by an 
order of magnitude for TeV photons. VERITAS should reach an angular resolution of 5 arcmin at 100 GeV [30] . As 
shown in the left panel of Fig. 2. the optimal beam size is 1.4 arcmin in the direction of M87. For a 6 arcmin angular 
resolution and a 50 GeV threshold, the significance drops down to 5 = 0.2. A 3-cr detection level translates therefore 
into the gamma-ray production cross section (crw)coiit. -^7 = 1-5 x 10~^^ cm^ s~^ for a 1 TeV neutralino. The heavy 
solid line in the upper-left panel of Fig. 3 corresponds to a sensitivity level of 

Mco„t.^7 > 1-5 X lO-'^cm^s-i { ^QQQ^y Y ■ (27) 

The region below that line will not be accessible, even with the next generation of Cerenkov telescopes. When the 
threshold increases to 100 and 250 GeV, the sensitivity respectively reaches down a level of 7.2 x 10"^'' and 2.8 x 10"^'^ 
cm^ s^^. Notice that even for a 250 GeV threshold, the number of background photons collected within 6 arcmin 
from the center of M87 amounts to ~ 4,600 particles. A 3-a signal corresponds to ^ 200 additional gamma-rays from 
that hot spot. For a 50 GeV threshold, \/Nhg ~ 1-4 x 10^ photons whereas Ng ~ 1, 100 photons. 

In Fig. 4, we present the results for the gamma-ray lines. We considered the two processes xx ~^ 77 a^nd xx ~^ ■^°7- 
The photon energy in the 27 process is clearly E^ = m-^, while in the Z'^'y reaction the photon energy is 

E^=m^-^^. (28) 
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FIG. 3. Annihilation rates in the continuum channels. Three thresholds are illustrated, Eth = 50, 100 and 250 GeV. 
Considering M87 as the source, the 3-a detection limits for exposures of 0.01 km^ yr are also presented. The regions below 
the higher heavy solid lines will not be accessible, even with the next generation of Cerenkov telescopes. The lower solid lines 
show the region of accessibility if annihilation rate is enhanced by a factor of 40 due to dumpiness. 



As in Fig. 3, we also present detection limits. The left panel features the two gamma-ray line whereas the right 
panel presents the Z^j process. The sensitivity limits are identical for both lines at high energy, i.e., for massive 
neutralinos. In the case of the Z'^j line, the photon energy becomes vanishingly small as the neutralino mass 
tends to mz/2. Because the signal becomes swamped inside a very strong low-energy gamma-ray background, the 
sensitivity drops completely and the detection limit of the right panel exhibits a sharp increase around mz/2. Below 
that threshold, the process is no longer kinematically allowed. There is an additional restriction arising from 
the energy threshold i?th of the Cerenkov telescope itself. Requiring that > Eth implies that the neutralino mass 
should exceed 
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FIG. 4. Annihilation rates in the monochromatic channels, (left): 2-y line, (right): Z^'y line. Both processes are present 
with their S-cr detection limits for a 10 km^ yr exposure towards M87 together with an energy resolution of AE/E = 0.02. 



Eth + \A 



th 



(29) 



For a 50 GeV threshold, this translates into > 77 GeV. We find that the lines are much more difficult to detect. 
The sensitivities presented here correspond to an exposure of 10 km^ yr and an energy resolution AE/E = 0.02, both 
of which are unreasonable with today's detectors. 



VI. DISCUSSION AND CONCLUSIONS 



As is clear from Fig. 3 and 4, the supersymmetric parameter space will mostly remain below the sensitivity level of 
the next generation of instruments. A few configurations are potentially detectable provided that the ACT threshold 
is decreased down to 50 GeV. If neutralinos are clumped inside the dark matter halo of M87, the situation considerably 
improves. As the gamma-ray production rate goes as the square of the density, neutralinos annihilate more efficiently 
as they condense. That overall effect may be described by the dumpiness factor C, which is defined as the increase of 
the global annihilation rate that results from a possible clumpy structure of the dark matter distribution. 



Neutralinos are cold dark matter species. As such, they exhibit density fluctuations 

P J 



that are related to the comoving wave vector k through [31] 
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Ak 



(l + a/s + /3fci-5 + 7fc2)2 



(30) 



(31) 



with a = 1.71 X /3 = 9 X l^-"" and j = f where / = (fJft.^)-!. Normalization to as = 0.8 gives A = 2.82 x 10^ Mpc'' 
when = 1 and h — 0.5. In a restricted wavelength range, it is approximated by a power law 



P{k) cx F 



(32) 



The power spectrum P{k) of density fluctuations behaves as fc on small scales, i.e., for structures typically lighter 
than 10® Mq. As regards a possible clumpy structure of the halo around M87, the relevant mass range extends from 
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Mi ^ 10^ M© up to Ms ~ 10"'^^ Mq. The corresponding spectral index n goes from —2.6 to —2.1. As shown below, 
structures smaller than M{ turn out to all have the same density. They contribute identically to the dumpiness factor 
C. There is no larger structure than the halo itself whose mass Mg reaches 10^'^ Mq in the inner 100 kpc. Because 
the comoving wave vector k scales as M~^/^, neutralino density fluctuations depend on both the scale M and the 
redshift z as 

^ oc (1 + z)-^ M-("+3)/6 . (33) 
p 

The redshift factor (1 + z)^^ is typical of the t^/^ growth of density fluctuations in a flat matter-dominated universe. 
Notice that small scale pcirturbations, for which n = —3, all become non-linear at the same time. Their subsequent 
collapse leads to virializcd structures whose densities have been enhanced by a factor of ~ 180 with respect the epoch 
of formation, when 5p/ p reached unity. Small scale dark matter clumps all have therefore the same density today. 
The formation redshift of larger structures behaves as 

(1 + ^f) oc M-("+3)/<' , (34) 

so that today, neutralino clumps with mass above ~ 10^ Mq have a density 

p{M) cx 180 (1 + z^f oc M-("+3)/2 . (35) 

The density p{Ms) of the largest possible clump should be comparable to the average dark matter density Pdm in the 
halo around M87. The distribution of clumps should follow the Press-Schechter's law 

dN _ Mo 

dM ~ M^ ' ^•^^> 

The normalization mass Mq obtains from the requirement that the clumps make up a fraction / of the halo. Disre- 
garding for the moment clumps with mass less than 10^ Mq, we get 

Some clumps arc actually destroyed through the tidal stripping resulting from both their mutual interactions and the 
action of the gravitational field of M87. Just like globular clusters orbiting the Milky Way, they evaporate so that a 
fraction / only of the initial population is expected to survive. Inside a clump with mass M, the annihilation rate of 
neutralinos is {av) {p{M) / m^} per unit volume. A net number {av) p{M) M/m^ of annihilations take place in the 
clump per unit time. We infer that the total annihilation rate of clumped neutralinos is obtained from the convolution 

Tciump = / ^p(M)M^rfM. (38) 

jMi 

Taking into account the mass-density relation (35) as well as the mass distribution (36) of the clumps, we readily 
infer the rate 

rau.p = ^P-Mo/^^ (-) - . (39) 

This may be compared to what a homogeneous distribution would yield 

Thorn = -^PdmM, . (40) 

ml 

The dumpiness factor C may be understood as the enhancement ratio rdump/rhom- Clumps span less space than 
if their matter was homogeneously distributed. In their interiors, neutralinos nevertheless annihilate much more 
efficiently. The net effect is the increase 



C = 



f 



n + 3j V In (Ms/Mi) 



\Mi ^ 



(41) 
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This expression gives C w 34/ for n = —2.1 and C w 4/ for n = —2.6. With a general power spectrum P{k), this 
expression reads 



C = (^^^] r (^] (^] '-^ (42) 



^ In (Ms/Mi) ; Jm, \MsJ \P{Ms)J M 
where we go from P{k) to P{M) through 



^=(ri^^Tj2-2 ) " (43) 



M G 



The CDM power spectrum leads to C « 13/. If we now assume that most of the clumps are small and that their 
mass does not exceed Mi = 10^ M©, we find 

Tciump = ■^p(Mi) f MdN . (44) 
ml J 

If that population of light clumps accounts for a fraction / of the dark matter halo around M87, the previous relation 
translates into 

:^ciump = -^p(Mi)/Ms . (45) 



The dumpiness factor becomes 



l,{M,) , (M. - 1"«>'^ 



in the case of a power-law spectrum or more generally 



\MJ \P{M, 

It reaches a value of C ^ 40/ for the CDM power spectrum. Varying the fraction / between 0.1 and 1, wc conclude 
that depending on the typical size of the clumps, the gamma-ray production rate may be enhanced by factors as large 
as 40. The lower solid lines in Fig. 3 show the sensitivity limits of ACTs assuming that C = 40. A Cerenkov telescope 
operating with a 50 GeV threshold would detect a neutralino-induced gamma-ray emission from the giant elliptical 
galaxy M87 for a part of the supersymmetric configurations outlined in the upper-left panel of Fig. 3. Even with the 
annihilation rate enhanced by a factor of 40, the gamma ray lines are out of reach. 
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